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Abstract 

The so-called semi-relativistic expansion of the weak charged current in powers of the initial nucleon 
momentum is performed to describe charge-changing, quasielastic neutrino reactions at 

intermediate energies. The quality of the expansion is tested by comparing with the relativistic 
Fermi gas model using several choices of kinematics of interest for ongoing neutrino oscillation 
experiments. The new current is then implemented in a continuum shell model together with 
relativistic kinematics to investigate the scaling properties of (e, e') and {y^,yL~) cross sections. 
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1 Introduction 

The importance of neutrino-induced reactions in nuclei has been stressed in connection with 
the neutrino oscillation experiments performed by the KARMEN and LSND collaborations 
[l]-[8]. In experiments of this type relatively low z/g or energies are involved (at most a 
few hundreds of MeV) and so the nuclear excitations involved can be described by standard 
non-relativistic models of the reaction including the relevant machinery (RPA correlations, 
large-basis shell models, A-hole excitations, hnal-state interactions, etc.) for this energy 
regime, where in particular giant resonances may play an important role [9]-[16]. 

However, when passing to the ongoing and next generations of neutrino experiments, 
MiniBooNE, K2K/T2K, MINOS, NOz/A, and MINERz/A [17]-[22], the neutrino beam ener¬ 
gies increase to the GeV level, and typically large energies and momenta are transferred to 
the nucleus. For these kinematics relativity is important and non-relativistic models of the 
reaction such as those listed above are bound to fail unless the relevant relativistic ingredients 
are included. 

First of all the use of relativistic kinematics is required, and must be implemented in the 
model. From quasielastic (QE) electron scattering studies we know that a good approxima¬ 
tion to the correct kinematics consists in the substitution 

A-^A(1 + A) (1) 

with A = u/2mM, where u is the energy transfer and uin the nucleon mass. This substitu¬ 
tion can easily be performed in all places in the calculation — except in the nucleon form 
factors, where the correct value of momentum transfer = [u, q) must be used. Second, 
a good approximation to the current matrix elements is required. The current is tradition¬ 
ally obtained from the fully relativistic one by some expansion procedure, usually through 
the Foldy-Wouthuysen transformation, which is valid for small momenta compared with the 
nucleon mass. Within this procedure relativistic effects were studied for neutrino energies 
up to 300 MeV in [23], showing that a good agreement is obtained between fully relativistic 
and non-relativistic Fermi gas calculations, when terms up to order are included 

in the latter. Similar results were obtained at larger energies, for the scattering of atmo¬ 
spheric neutrinos from oxygen [24], where, however, it must be noted that, due to the steep 
decrease of the atmospheric neutrino spectrum at large energies [25], the most important 
contributions to the process correspond to relatively small energy transfers. Obviously the 
Foldy-Wouthuysen expansion is not applicable for the high values of the momentum transfer 
q which are of interest for the new experiments, i.e., for values around 1 GeV/c. In this 
case a different expansion procedure must be performed in which q and u can be arbitrarily 
large. 
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The first goal of this paper is thus to develop an approximation to the nuclear charged 
current (CC) that accounts for specihc relativistic effects which are relevant to intermediate 
energy quasielastic neutrino reactions. The CC is obtained here through an expansion that 
only requires the momentum of the initial nucleon to be small, while it treats exactly the 
dependence on (cu, q). Moreover the resulting semi-relativistic (SR) current is simple enough 
to be easily implemented in already existing non-relativistic models of and (l7;,/+) 

reactions, where I = e or fi. It is presented here as an extension of the electromagnetic and 
weak-neutral current expansion originally derived in [26] and that in recent years has been 
widely tested and applied in several collaborations to describe a wide variety of inclusive 
and exclusive electron scattering observables for intermediate energies and excitations in 
the vicinity of the QE peak [27]-[37]. Extensions of the SR expansion to meson-exchange 
currents have also been developed [38, 39], and a detailed description of their application to 
two-body currents can be found in a recent review article [40]. 

In this paper we apply the SR model to neutrino-induced reactions for the hrst time. 
Apart from the different isospin dependence, the SR expansion coincides essentially with 
that of the weak neutral current performed in [26]. The only difference is that in the present 
reaction we have to include the time component of the axial-vector current. That component 
was not considered in [26] because it does not contribute signihcantly to parity-violating 
electron scattering. 

We check the quality of the SR expansion in the context of the Relativistic Fermi Gas 
(REG) model. This is a very convenient model for our purposes, since it is fully relativistic 
and simple enough to be solved exactly. Also this model is capable of getting the basic 
size and shape of the QE (e, e') cross section. Thus we shall show that, starting from the 
non-relativistic Fermi gas, performing the replacement in Eq. (1) in the kinematics and im¬ 
plementing the new SR charged current, we reproduce basically the REG results for (z/^, ^~) 
reactions. This “relativizing ” procedure can easily be extrapolated to more sophisticated 
hnite nuclei models of the reaction. The result will at least reproduce appropriately the 
allowed kinematical region and the relevant relativistic content of the current operator. Of 
course features related to the relativistic aspects of the dynamics cannot be accounted for by 
our procedure. In context we note that studies of relativistic nuclear dynamics in charged 
and neutral-current neutrino-nucleus QE scattering have already been presented in some 
previous work [41, 42, 43, 44, 45]. In these studies a basic focus was to analyze the effects 
introduced by various descriptions of the hnal nucleon relativistic states upon the integrated 
cross sections. 

We illustrate the relativizing procedure by applying it to the continuum shell model 
(GSM), i.e., nucleons in a mean held taken here to be a Woods-Saxon potential. We use the 
same (real) potential for initial and hnal states in order to maintain orthogonality between 
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nuclear states. In this way we extend the SR shell model of [26] to neutrino reactions at 
the QE peak. Since the use of relativistic kinematics is potentially equivalent to solving a 
Klein-Gordon equation, the present SR continuum shell model includes some aspects of a 
relativistic mean held. 

With such a model we are in a position to fulhll the second goal of this paper, which is to 
investigate superscaling properties of both (e, e') and inclusive cross sections at the 

QE peak for intermediate energies, namely the degree to which one hnds that the reduced 
cross sections are independent of the momentum transfer (scaling of the hrst kind) or the 
nuclear species (scaling of the second kind) or both (superscaling). Exhaustive analyses of 
the (e, e') world data and explorations of various aspects of their scaling properties have been 
performed in [46]-[52]. In particular, in recent work [53] the approach has been extended 
to the A peak, allowing one to construct a semi-empirical model based on scaling which is 
very successful in describing the experimental (e, e') cross section up to the A peak for high 
energies. This model allowed us in [53] to generate predictions for the (z/^,/i“) cross section, 
under the reasonable hypothesis that it presents, for high energies, scaling properties similar 
to those of the (e, e') cross section. This hypothesis is true by construction in the REG and 
can be also demonstrated at least for the conventional Plane Wave Impulse Approximation 
(PWIA) in the high q limit However, when distortion of the ejected nucleon is present, a 
general proof of scaling cannot be provided, while for lower energies it is clear that the axial- 
vector and vector matrix elements are renormalized differently in the nuclear medium due to 
RPA correlations, and this can obviously modify the scaling properties of the neutrino cross 
section. For higher energies one expects that nuclear effects such as from RPA correlations are 
less important than at low energies. In this work we use the SR shell model to investigate the 
degree of violation of the scaling hypothesis for the QE peak region within the Distorted Wave 
Impulse Approximation (DWIA), where the distortion of the ejected nucleon is described 
with a real potential. We perform this study in two steps. First we focus on the (e, e') 
cross section and study the scaling properties of the separate response functions. Once the 
superscaling has been verihed for electromagnetic processes, we are able to reconstruct the 
(z/^,/i“) cross section from the (e, e') one using the scaling hypothesis and compare it with 
the one computed directly using the SR shell model. Thus, at least within the context of the 
SR shell model discussed in the present work, following this procedure we shall be able to 
check the consistency of our approach and quantify the degree to which scale-breaking effects 
are expected to enter. As we shall see below, there appears to be very little impact from 
this source of scale-breaking on the scaling approach used in [53] to predict neutrino-induced 
cross sections. 

^However, this is broken to some degree by the Relativistic PWIA (RPWIA) where factorization no longer 
obtains. 
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The paper is organized in the following way. In Sec. 2 we begin with a brief review of the 
general formalism for neutrino scattering, and present the expansion of the CC. We particu¬ 
larize the formalism for the shell model, and introduce the general multipole expansion of the 
responses, with some details on the derivation of the Coulomb multipoles of the axial-vector 
current placed in Appendix A. We also provide the expressions for the factorized PWIA in 
Appendix B. In Sect. 3 we present results for the reaction for several choices of 

kinematics and several nuclei of interest. We hrst perform an analysis of the quality of the 
relativizing procedure by comparison with the RFC for the relevant kinematics. The quality 
of the various components of the SR current can be checked separately by examining the 
individual response functions that contribute to the process. We then focus on the SR shell 
model and perform the scaling analysis of the longitudinal and transverse electromagnetic 
responses as functions of the momentum transfer for various nuclei. We apply the scaling 
hypothesis to reconstruct the cross section starting from the electromagnetic scal¬ 

ing function and compare with the shell model results. Finally, in Sec. 4 we present our 
conclusions. 


2 Formalism 

In this section we briefly present the basic formalism for neutrino-induced reactions. Some 
of the previous approaches to the general formalism for these reactions can be found in [41], 
[54]-[57]. 


2.1 Charge-changing neutrino cross section 


Here we focus on the particular case of while the cases of anti-neutrinos or of 

other lepton species can be easily obtained with obvious changes. The four-momenta of 
the incident neutrino and detected muon are = (e, k) and k'^ = (e',k'), respectively. 
The four-momentum transfer is = k^ — k'^ = (cc;,q). We use a coordinate system with 
the 2 :-axis pointing along q and the x-axis along the transverse component of the incident 
momentum, i.e., k-*“ = k — ^(k'q)q. We follow the formalism of [53] where the cross section 
is written as 


da 

dkl'de' 


aoJ^l 




( 2 ) 


with 


cos^ Oc, , / 2 ^ 

an = -r —k e cos - 

“ 27r2 2 


( 3 ) 
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Here G = 1.166 x 10 MeV ^ is the Fermi constant, 6c is the Cabibbo angle, cos^c = 
0.975, and the angle 9 is dehned as 


2 ^ 

tan - = 




2 (e + e'y — 

with Q"^ = <{). The nuclear structure information is contained in JFJ, dehned by 


( 4 ) 


= VccRcc + ‘^VclRcl + VllRll + VtRt + 2Vt>Rt' , (5) 


where the kinematical factors Vk coming from the leptonic tensor are dehned by 


VcL 

Vll 

Vt 

Vt' 


1 — (5^ tan^ - 
2 


uj (5^ 


■ tan 


9 


G 2U 2 \ r 2 +„„2 ^ 


= tan" 


q p' 

+ (1 + ft. + ) (5^ tan^ 

q^ V (IP ) 

9 p 5 “^ (oj 
2 + 2 “ 7 
J 




2£_ 

d' V <1 2 


In Eqs. (6-10), following [53], we have dehned 

5 = ■ 

P = ■ 


m 


P = 


m 

q2 

q 

6 T 6^ 


( 6 ) 

(7) 

( 8 ) 
( 9 ) 

( 10 ) 


( 11 ) 

( 12 ) 

(13) 


Note that the only dependence on the muon mass m' is contained in the <5 coefficient. 
Finally the weak response functions are given by 


Rcc = 


(14) 

Rcl = 

_1 (w03 ^ ^30^ 

(15) 

Rll = 


(16) 

Rt = 


(17) 

Rt' = 


(18) 
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in terms of the inclusive hadronic tensor [54]: 


fi 

In Eq. (19), J^{Q) is the hadronic CC current operator, to be specihed below, and a sum 
over hnal states and an average over initial spin is assumed. 


2.2 Semi-relativistic charge-changing current 


We begin with the basic relativistic charged weak current of the nucleon, = jy — j^. In 
this work we employ only the Standard Model of electroweak interactions at tree level and 
thus, for example, do not include radiative corrections or contributions from second-class 
currents (see [54]). We use the conventions of [58]. The vector and axial-vector currents are 
given by 


Jv(p,p) = w(pO 
ja(p,p) = w(pO 


2 FiV+*— 

rriN 


«(P) 


+ Gi 




2mM\ 


7^m(p) , 


( 20 ) 

( 21 ) 


where for the isovector nucleon form factors EY 2 = {F 12 ~ use the Galster pa¬ 

rameterization [59] and m(p) is the free Dirac spinor of the nucleon. The axial-vector and 
pseudoscalar form factors are parameterized as 


Gp 


9A 

1 - QVAQ 

A 


( 22 ) 

(23) 


with qa = 1.26, Ma = 1032 MeV. 

The SR approximation to this current is then obtained by inserting the appropriate free 
spinors m(p), m(pO) and 7 matrices, and performing an expansion in powers of rj = p/ruN 
to first-order. The procedure was developed in [26, 38, 40] for the electromagnetic and 
transverse neutral, vector and axial-vector currents, and is based on the fact that inside the 
nucleus \r]\ is a good expansion parameter, the characteristic dimensionless nuclear scale be¬ 
ing Tjp = kp/mpf ~ 1/4. We exploit QE kinematics, while further reasonable simplihcations 
are eventually needed in order to arrive at simple expressions that are easily implementable 
in traditional non-relativistic calculations. For the vector current we use the following SR 
approximation; 


Jv = ^o + ^^o(«; X 77)-cr 

Jy = 6^7^ + X K , 
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(24) 

(25) 



where 


, ^0- 


H _ V ^ f _ r\/~iV \ 

SI — ) SI — 

K K 

and use has been made of the dimensionless variables n = q/ 2 m 7 v and r = — \^. 

From vector current conservation, the longitudinal component is given by Jy = ^Jy- 
Note that in Eqs. (24,25) the various terms making up the current are similar to the ones 
that can be found in traditional non-relativistic expansions commonly used for the charged 
current (see for instance [16]), except for the k and r-dependent factors, that provide 
the required relativistic behavior. In Jy we include the first-order (0(7^)) contribution. 
This spin-orbit term is proportional to the operator {n x rf) ■ cr and is of some importance 
for high q values. On the other hand, the transverse component Jy is the sum of the usual 
magnetization (crx k) piece plus a first-order term, the convection term which is proportional 
to 77 -*-, that gives in general a very small ^ contribution to the cross section for high q [26, 28]. 

In the case of the axial-vector sector, only the transverse component of the weak neutral 
current was expanded in [26]. We use the following version of the corresponding SR current 

Ji = c: = . (28) 


(26) 

(27) 


Note that we have neglected the terms of order rj since, as we shall show in the next section 
(and also demonstrated in [26] in the context of parity-violating electron scattering), they are 
small while adding unnecessary complications to the shell model calculation. They can be 
safely neglected for our purposes. Remarkably the factor yjl + r in Eq. (28) already accounts 
for the most part for relativistic effects in this current (see [26] for the full expansion of this 
current to hrst-order in rj). 

We are left with the 0 and components of the axial-vector current. Their SR expressions 
are presented here for the hrst time. Using a notation reminiscent of that used in Appendix A 
of [26], they are written to first-order in rj as 

Ja = Co«^ • o- + ■ cr (29) 

Ja = (s'* ■ o- + ■ o- , (30) 


where 





C" 

80 

C" 

83 




Ga- 


Ga- 




+ 1 ) 


Gh 


K 


K + + 1 ) 


Gh 


(31) 

(32) 


^Note that 77 


p/niN is essentially the velocity v of the initial struck nucleon in units of c. 
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and we have introduced the following combination of axial-vector and pseudoscalar form 
factors 

G'^ = Ga- tGp . (33) 

For these components of the axial-vector current, Eqs. (29,30), we have performed the ex¬ 
pansion to hrst-order in rj. The first-order axial-convective term is proportional to cr ■ r]^, 
and, as for the spin-orbit and convection terms, only the perpendicular velocity appears. 
As we shall show below, for the kinematics of interest in this work the form factor which 
drives the zeroth-order terms turns out to be small due to cancellations in Eq. (33) between 
Ga and rGp. In such cases the 0{r]) term is dominant in these current components. 


2.3 The continuum shell model 


In this work we restrict our attention to the closed-shell ^ nuclei and ^°Ca, which we 

describe in a Continuum Shell Model (CSM). Thus the present model does not include nuclear 
correlations. Such effects are important for low energy, in particular in the axial responses, 
and have been estimated for instance in the RPA approach of Ref. [15], but are expected to 
be smaller at the GeV energies considered in this paper. The initial state |i) appearing in the 
hadronic tensor, Eq. (19), is described as a Slater determinant representing the uncorrelated 
nuclear core with all shells occupied. Since we are working in the impulse approximation, 
the hnal states are particle-hole excitations coupled to total angular momentum J, namely 
I/) = \{ph~^)J). The single hole, \h) = \ehlhjh), and particle, \p) = \eplpjp), wave functions 
are obtained by solving the Schrodinger equation with a Woods-Saxon potential 


R(r) 


~^/(a -Rq; no) + 


Vis df{r,Ro,ao) 
m^r dr 


I rr + Vc{r) , 


(34) 


where 

f{r,R,a)= (35) 

and Vc(r) is the Coulomb potential of a charged sphere of charge Z — 1 and radius Rq (it is 
equal to zero for neutrons). The parameters of the potential are htted to the experimental 
energies of the valence shells and are given in Table 1. In the shell model the energy transfer 
is computed as the difference between the (non-relativistic) single-particle energies of particle 
and hole u = ep — eh- The relativistic kinematics are taken into account by the substitution 


ep^ ep{l + ep/2mN) 


( 36 ) 


as the eigenvalue of the Schrodinger equation for the particle; c.f., Eq. (1). 


^By shell here we mean a sub-shell with quantum numbers (nlj). 



VEs ro ao 


Kf Vis VI 


12C 

62.0 

3.20 

60.00 

3.15 

1.25 

0.57 

160 

52.5 

7.00 

52.50 

6.54 

1.27 

0.53 

4°Ca 

57.5 

11.11 

55.00 

8.50 

1.20 

0.53 


Table 1: Woods-Saxon potential parameters for protons (p) and neutrons (n). The units are MeV 
for Vi, and fm for oq and tq. The reduced radius parameter tq is defined by Rq = . 


Since the nuclear states have good angular momentum, it is convenient for the shell model 
calculation to perform a multipole expansion of the components of the current operator in 
terms of the usual Coulomb {C), longitudinal (L), transverse electric {E) and transverse 
magnetic (M) operators [26, 54, 60], dehned by 


C'jo(g) = J d^rjj{qr)Yjo{r)Jo{r) 

(37) 

Ljo{q) = - f d^rV[jj{qr)Yjo{r)]-J{r) 
q J 

(38) 

Ejmiq) = - f X [jjiq'r)'^ JJm{i)] ■ J(r) 

q J 

(39) 

Mjmiq) = y dVjj(gr)Yjj^(f )■ J(r) , 

(40) 

where jj{qr) is a spherical Bessel Function and Yjjm(f) is a vector spherical harmonic. 

The nuclear response functions are then written as 

Rcc = 4^E(|Clt + |C^r) 

a 

(41) 

Rcl = 27rY.{cl‘Ll+ClLl' + Ci-Li + CtLf) 

a 

(42) 

rll = 4^i:(iiri"+ii^i0 

a 

(43) 

Rt = i7rY.{\El? + + \Et? 

a 

(44) 

Rt' = 27r ^ (El* Ml + El Ml* + El* Ml + El Ml*) 

a. 

(45) 

where in the sums Cp = e/j + cu is modified according to Eq. (36), and we use the index a 
to label the quantum numbers a = {h,lp,jp, J). Moreover, the C, L, E and M multipoles 
of the vector and axial-vector currents are dehned by the reduced matrix elements of the 
corresponding operators 

Cl + ^Cl = {f\\Cj{qm 

(46) 
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Ll + iLt 

= {fwUm 

(47) 

El + lEl 

= (fWMgW) 

(48) 

. 

1 

= {fWMAqm. 

(49) 


These reduced matrix elements are given in [26] for the vector and axial-vector transverse 
operators, and in [60] for the leading-order longitudinal axial-vector component. In the 
present paper we add the hrst-order convective term of the axial-vector current ■ a 
appearing in Eqs. (29,30). The and multipoles of this new term are presented in 
Appendix A. 

The sums over a in Eqs. (41-45) are inhnite and have to be truncated in the calculation 
once convergence is reached. In our approach the number of multipoles is determined by 
the maximum of the total angular momentum, Jmax- We £x this quantity by computing the 
response functions setting the potential in the final state to zero and comparing with the 
factorized PWIA (see Appendix B for details). The number Jmax increases with q and with 
the number of nucleons, A. In the next section up to 41 multipoles have to be summed up 
for the case where g = 1.5 GeV/c. 


3 Results 

3.1 Test of the SR approach 

The quality of the SR expansion of the charged current is illustrated in Figs. 1-3. There we 
show the separate vector, axial-vector and total response functions of for three values 
of the momentum transfer, q = 0.5, 1, and 1.5 GeV/c. The Fermi momentum is chosen 
to be kp = 220 MeV/c. We show the separate response functions instead of cross sections 
because in these functions we can appreciate better the contribution of the various current 
components and the behavior of the SR approximation as a function of the momentum and 
energy transfers. This can be seen in the hgures by comparing the solid lines, represent¬ 
ing the REG result (exact relativistic answer) with the dashed lines, corresponding to the 
Semi-Relativistic Fermi Gas (SRFG). This last model is obtained by implementing the SR 
expansion of the current in a non-relativistic Fermi gas, although with relativistic kinematics. 
The accord of the two models is almost perfect for the three values of q considered. In Fig. 1 
we only show and since and related to by current conservation. 

In Fig. 2 we show instead the three response functions, since the axial-vector current is not 
conserved. Finally, in Fig. 3 we show the sums of Figs. 1 and 2 and also the interference 
response Rp'. 

In Figs. 1-3 we also show with dotted lines the GSM responses, computed using the 
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Figure 1: Vector response functions of for three values of the momentum transfer. Solid 
lines: RFG with kp = 220 MeV/c. Dashed lines: SRFG. Dotted: GSM. Dot-dashed: PWIA. 


relativizing procedure. For comparison we also show the PWIA results obtained by setting 
the potential in the hnal state to zero (or, equivalently, by integration of the factorized ex¬ 
clusive responses, as shown in Appendix B). The GSM response functions are quantitatively 
similar in magnitude to the RFG and are centered approximately around the same value 
of OJ. The major discrepancy between the two models is found for the axial-vector Rcc) 
and Rf^j^ response functions, where the GSM responses are slightly larger in magnitude 
than the RFG ones. These responses are in general small compared with the corresponding 
vector responses. This is a consequence of the suppression of the axial-vector current in the 
longitudinal channel due to small value of the form factor dehned in Eq. (33), for these 
kinematics. In fact, from the dehnition in Eq. (23) of the pseudoscalar form factor one has 

At the intermediate energies of interest, \Q‘^\ 3> and hence the factor inside the parenthe¬ 
sis in Eq. (50) is also small, of order 0{m‘l/Q‘^). In this situation the hrst-order, axial-vector 
convective term of the current, which is proportional to ry*- -cr (see Eqs. (29,30)), is dominant 
over the zeroth-order contribution, and the corresponding response functions are in general 
small compared with the vector ones. An example is shown in Fig. 4 where we display the 
separate contributions of the zeroth- and hrst-order terms to the axial-vector response Rqq 
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Figure 2: As for Fig. 1, but now for the axial-vector responses 






[GeV-1] [GeV i] Rll [GeV i] [GeV~^] Rcc [GeV-^ 



0 100 200 300 400 200 400 600 800 600 800 1000 1200 

Lo [MeV] to [MeV] to [MeV] 


Figure 3: As for Fig. 1, but now for the total (vector plus axial-vector) responses 
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of for q = 500 MeV/c in the CSM. For higher values of q the zeroth-order contribution 
is much smaller than the others. 



Figure 4: Axial-vector response of in the CSM. Solid: contribution of the zeroth- 
order term in the SR expansion. Dashed: hrst-order contribution. 


From inspection of the responses Rcc, Rt and Rt' we observe that the PWIA results 
are clearly shifted to the right of the CSM by roughly the averaged depth of the potential 
~35 MeV. This shift is present also in the separate vector and axial-vector responses. In 
the case of the CL and LL the shift is larger due to the energy factors used to compute 
the L component of the current. The origin of the shift in PWIA is related to the different 
treatment of the nuclear hamiltonian in the initial and hnal states [27]: while the energy of 
the initial bound neutron is the sum of kinetic plus potential, = th + Vh-, the exiting proton 
is a plane wave and has only kinetic energy Cp = tp. Accordingly, for hxed u there is an 
imbalance between potential energies in the initial and hnal states, yielding uj = tp — th—Vh- 
However in the CSM we use the same potential for the initial and hnal state, which now 
has ep = tp + Vp. Hence u = tp — R + Vp — Vh- Thus in the CSM the potential energies of 
particle and hole partially cancel out and this explains why the position of the peak is close 
to the RFC, where only kinetic energies enter. Of course the cancellation is not perfect, and 
a slight shift to the right of the RFC is also observed in the CSM. This small shift of the 
CSM is the behavior expected from previous theoretical studies [61]. 

Another test of the SR approximation is illustrated by the results shown in Fig. 5 for the 
differential cross section ^^C(z/p,/i“). Therein we show examples for two incident neutrino 
energies, e = 1 and 1.5 GeV, and for two scattering angles 9 = 45°, 135°, as a function of 
the exiting muon energy e'; so in this case we are testing different ranges of q and u, which 
can be high or low depending on the kinematics. We see again that the RFC and SRFG 
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Figure 5: Differential cross section of the reaction for incident neutrino energies 

e = 1 and 1.5 GeV and for two scattering angles. Solid lines: RFG with kp = 215 MeV/c. 
Dashed lines: SRFG. Dotted: GSM. 


predictions are almost equal in all the cases, while the GSM also gives similar results, with 
the exception of the characteristic tails and slight shift, that now is to the left of the RFG 
since ui decreases with e'. 

Summarizing this subsection, using the Fermi gas as “testing arena” for the approximated 
GG current, we have found that the accord between the results obtained using the SR 
expansion and the exact relativistic result is almost perfect for the intermediate and high 
values of q and oj considered. The GSM relativized using our procedure gives rise to cross 
sections which are within the allowed kinematical region and of a magnitude similar to those 
of the RFG. This behavior of the GSM is similar to what was found for (e, e') reactions 
in [26]. We should underscore the fact that the GSM is used here just as illustration of 
how one can use our relativizing procedure in a more complex model than the Fermi gas, 
and of the results that are to be expected from it. In particular, especially when discussing 
lower energies where non-relativistic approaches are valid, more elaborated models exist that 
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describe the (e, e') experimental data rather well. For example, it is well known from other 
work that hnal-state interactions can signihcantly modify the bare CSM or RPA responses, 
mainly through medium renormalization of the particle-hole excitations. This mechanism 
can be approximately taken into account by using an effective nucleon mass that produces 
a shift and by a folding of the bare responses, producing an asymmetric broadening which 
improves the agreement with the experimental data [62, 63, 64]. 

3.2 Scaling 

In this subsection we work within the CSM. We hrst focus on the inclusive electron scattering 
reaction, and investigate the scaling properties of the electromagnetic responses for various 
choices of kinematics for and ^°Ca. The (e, e') cross section is given by 

dcT 

~ ^Mott [VlRl + VtRt) ■ (51) 

The same SR expansion for the vector current, Eqs. (24,25), is employed here for the elec¬ 
tromagnetic sector. 

We compute the CSM scaling functions as 

h = ^ (52) 

h ^ ^ (53) 

with 

Gk = Ao {ZUl + NUl) K = L,T, (54) 

where Aq is given in Eq. (94) and the electromagnetic single-nucleon functions for protons, 
and neutrons, U^, are dehned similarly to the ones given in Appendix B for the vector 
current, Uqq and . 

The scaling behavior can be studied by plotting these functions against the scaling vari¬ 
able [46] 



for different kinematics and for different nuclei. Here = y 1 -|- k‘p/m% — 1. 

A study of the behavior of Jl, computed for various values of the momentum transfer, 
is summarized in Eig. 6. We plot together /l for q = 0.5, 0.7, 1, 1.3, and 1.5 GeV/c. All of 
the curves approximately collapse into one. Small violations of the scaling are seen at low 
Ip coming from the low-energy potential resonances for q = 0.5 GeV/c, which disappear for 
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tp 

Figure 6: Scaling of the first kind in the CSM with scaling functions obtained from the 
longitudinal electromagnetic responses in the CSM. In each panel we include q = 0.5, 0.7, 1, 
1.3 and 1.5 GeV/c. 
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Figure 7: As for Fig. 6, but now including in each panel the nuclei and ^“^Ca. 
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Figure 8: Scaling function /x, obtained from the longitudinal electromagnetic response (solid 
lines), compared with the scaling function /x obtained from the transverse response (dashed 
lines), for several values of q. 
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higher g-values. Thus scaling of the hrst kind is approximately achieved in the CSM. Note 
that scaling also holds for \'ip\ > 1, the region where the RFG responses are zero. 

Scaling of the second kind, i.e., independence of the nuclear species for hxed g, is illus¬ 
trated in Fig. 7. We plot together for the three nuclei studied. The htted Fermi momenta 
are 220, 215 and 240 MeV/c for and ^°Ca, respectively. The collapse of the three 

curves into one is clearly seen, with small deviations at the region of the maximum. Again 
the exception is found for q = 0.5 GeV/c in the resonance region. Since both kinds of scaling 
are found, we conclude that superscaling occurs within our model. 

Finally, in Fig. 8 we show what has been called scaling of the zeroth kind [50], i.e., the 
longitudinal and transverse scaling functions, and Jt, also collapse into one universal 
function /. Experimentally, deviations from scaling in the region of the QE peak mainly 
occur in the transverse response, and are related to contributions beyond the impulse ap¬ 
proximation, in particular to meson-exchange currents [49, 50]. Since these contributions are 
not included in our GSM, the L and T responses scale in the same way, as shown in Fig. 8. 

In the present work we do not compare with the experimental (e, e') scaling function 
because, as stated above, the GSM is still lacking some ingredients (for instance, medium 
modihcations to the p-h propagator, inclusion of energy-dependent potentials, etc.), without 
which one should not expect to obtain excellent agreement with the experimental data. The 
focus of the present study has been limited to showing that superscaling occurs at the level 
of relativized GSM. 

Let us now turn to neutrino reactions. In [53] a semi-empirical model, based on the 
superscaling property of the L responses of ^^G and for high momentum transfer, was 
proposed to predict the neutrino inclusive cross sections up to the A peak. The latter 
were calculated starting from the RFG expression for the cross sections and substituting the 
“theoretical” RFG superscaling functions with phenomenological ones (one for the QE and 
a different one for the A peak), derived from hts of electron scattering data. Therefore this 
reconstruction, besides assuming the validity of superscaling in electron scattering, relies on 
the hypothesis that it also holds for the neutrino inclusive cross section for the high energies 
involved, so that electron scattering results can safely be used as input in neutrino scattering 
calculations. This assumption obviously cannot be tested using the RFG model, since, as 
stated in the introduction, it is true by construction. 

In the present work we have at hand a model, the GSM, for both the QE (e, e') and 
reactions that has been relativized and so should be able to handle modeling at 
high energies. Moreover, superscaling of the electromagnetic responses is well satished by 
the GSM. Hence we can adopt the same approach as in [53] treating the GSM electromagnetic 
scaling function as a pseudo-phenomenological one and using it to compute neutrino cross 
sections. Upon comparing cross sections obtained this way with those obtained directly 
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using the CSM we can then gain some insight into the degree to which the shell model 
incorporates effects which are scale-breaking. While the level of scaling-violation is expected 
to be quite small, since we have already seen excellent superscaling in the hgures discussed 
above, the specihc differences in the roles played by the various current operators involved in 
electromagnetic and CC weak processes might lead to different sensitivities to scale-breaking 
effects. Accordingly, it is useful to compare the cross sections obtained within the scaling 
approach with those computed directly using the model. In context it should be noted 
that this does not imply in general that the scaling approach has been shown to be robust, 
since here the comparisons are being made entirely within limited models which clearly lack 
elements that may play some role in the responses and may be scale-breaking — we know 
of at least one such ingredient, namely, meson-exchange currents which have been shown to 
provide scale-breaking corrections that enter typically at the 10% level in the overall cross 
sections. Nevertheless, the stndy presented below indicates that the scale-breaking effects 
incorporated in the present models are in fact qnite small. 

We show typical results of this study in Figs. 9 and 10 for and respectively. The 
cross section for (z/^, fi~) is shown for two incident energies and for two scattering angles. The 
dashed lines have been compnted from the scaling function fi obtained from the analysis of 
(e, e') CSM predictions. To be precise, we use Eq. (93) to compute the neutrino responses, 
substituting by fi. The solid lines correspond to the exact CSM result. The differences 
found between the two approaches are quite small, at most ~ 3%. These small differences are 
produced by the slightly different scaling behavior of the axial-vector and vector responses, 
that, however, has little effect on the total cross section. Due to the simplicity of the CSM, 
these results cannot be taken as a dehnite proof of snperscaling, bnt certainly represent an 
important step forward in establishing the validity of the approach presented in [53]. 

4 Conclusions 

One of the goals of the present stndy has been to explore the degree to which scaling and 
superscaling behaviors are reached for relatively high-energy semi-leptonic inclusive reactions 
with nnclei at excitation energies in the vicinity of the qnasielastic peak. The focus has been 
placed on comparisons between the Fermi gas model (both the fully relativistic Fermi gas 
and a semi-relativized version of it) and a relativized continnnm shell model. 

Since we are interested in energies of several GeV, relativity is clearly a required ingre¬ 
dient. In the present work we have presented a clear and direct procedure that allows one 
to incorporate some classes of relativistic effects and thereby to relativize otherwise non- 
relativistic models of the reaction, such as the shell model employed here. We use the two 
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e = 1 GeV 


e = 1.5 GeV 



Figure 9: Differential cross section of the reaction for neutrino incident energies 

e = 1 and 1.5 GeV and for two scattering angles. Solid lines: GSM. Dashed lines: recon¬ 
structed from the electromagnetic scaling function /l ('?/’) computed at the same kinematics. 


versions of the Fermi gas model to motivate the procedures followed. Two steps are involved: 
hrst, one relativizes the kinematics in the reaction, and second, expansions are made for the 
single-nucleon currents in the problem. In the present work we have extended our previous 
treatments to include the full charge-changing weak interaction current. Importantly, be¬ 
cause the expansions are made to hrst-order in p/rriN, but not at all in q/ruN or u/rriN, they 
differ from the traditional FW expansions and for QE scattering should be more robust at 
high energies where expansions in q/mj^ and u/ttim clearly fail. 

Using the two versions of the Fermi gas model and the relativized continuum shell model, 
together with some results from the plane-wave impulse approximation, to represent both 
electromagnetic, (e,e'), and GG weak, processes, we have proceeded to quantify 

the level of scaling violation. This allows us to evaluate the corresponding uncertainty in 
the predictions one makes for the neutrino cross section using the scaling approach with 
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Figure 10: As for Fig. 9, but now for 

input from analyses of electron scattering data. In the present paper we have shown that for 
intermediate to high energies the description of the hnal state through a real mean held seems 
not to change appreciably the scaling and superscaling properties of the electromagnetic 
response functions at the quasielastic peak, and that the scaling function extracted from 
these is basically the same as for the reaction. 

Although the present analysis performed within the relativized shell model does not 
constitute a proof of the scaling hypothesis in the general case, scaling studies such as those 
presented here allow us to gain additional insight and, in particular, to rule out (or not) 
some reaction mechanisms as possible causes of scaling violations. 

In summary, the relativizing procedure followed in the present work embodies some of 
the ingredients that are almost certainly required for a sucessful description of such inclusive 
processes in the quaselastic regime — models which do not take at least these ingredients 
into account are very likely to fail in the several GeV energy region of interest here. This 
does not mean, however, that additional dynamical features are not needed before a full 
understanding of these high-energy responses will be attained. In particular, other studies 
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being pursued in parallel by us hold some promise for reaching a good understanding of 
the phenomenologically-derived scaling function used recently to make predictions for CC 
neutrino reactions in the GeV region. Some of these results will be presented in the near 
future. 
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A Multipoles of the convective axial-vector operator 

Here we compute the Goulomb multipoles of the convective axial-vector charge operator, 
which is the hrst-order term in the expansion of (Eq. (29)). The corresponding coordinate- 
space operator is given by 

Pc(q) = -cr , (56) 

where r is the coordinate of the nucleon over which the operator acts. To obtain the multipole 
operators we proceed by writing the above operator in the form: 

pc(q) = (v - 

= O ■ ■ (58) 

Performing in this equation the multipole expansion of the plane wave for q along the ; 2 -axis 

= v^^i^[J]jj(gr)ljo(f), (59) 

j 

where [J] = ^/2J + 1, we obtain 

pc(q) = , (60) 

J 

where the Goulomb multipole operators are 

Cjo{q) = -[ 1 + jM'r)YjQ{i) . (61) 
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Later on we will also make the substitution p ——iV, although in this equation V operates 
only on jj{qr)Yjo{r). 

The procedure now is to use Racah algebra and the general properties of spherical har¬ 
monics and spherical Bessel functions to write the above operator in explicit spherical-tensor 
form. The following expressions can be obtained: 


C'jo 

= CoF. 

'-'JO 

i 1 
m [J] 

'-'JO 

i , 

J + s, J" 

= J' + s' 


_i_ 

'JO “1“ '-"JO 


'JLJ 


^ ,S. jlTi • 


(62) 

(63) 

(64) 


function, which is a linear combination of the Bessel functions for J±l), and we have dehned 
the auxiliary coupled operator 


= jj"(qr) [a ® |lj"(f) 0 V] 


(65) 


Next we proceed to compute the reduced matrix elements of this R-operator between 
shell-model particle and hole states \p) = \^lpjp), and \h) = \^lhjh)- Using standard Racah 
algebra [65, 66] we obtain 


{p\\Uj>'j'j\\h) 




If. 


MjMWV] 


2 Ip jp 
2 jh 
1 J' J 


X 


X! [Lh]sh{lh + 

s;i=±l 


J" 1 J' \l Ip J" Ln \ 

ih Ip 0 0 J 


X 



drr‘^Rp{r)jjff{qr) 


( d lh + Ss^-i\ 


Rh{r) , 


( 66 ) 


where Rp{r) and Rh{r) are the radial wave functions, Lh = Ih + Sh, and Sh = ±1- We can 
now use this result in Eqs. (63,64). It is convenient to use the following identities involving 
products of six-j and three-j coefficients, for J" = J' + s' and Lh = lh + Sh, with s', Sh = ±1: 


J" 1 J' \l Ip J" Lh \ 
lh Ip Lh]\^ 0 0 J 


p+ 

TIMThl 


[{lh + <^s;i,-l) {J' + <^*',- 1 )] ^ 


lh 

1 


J' Ip \ 

-1 0) 


— ShS' [{lh + <^Sh,l) {J' + <^s',l)] ^ 


lh J' /p \ 1 

0 0 0 yj 


(67) 
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f J 1 

J 1 

i 

J 

Lh \ 

^h Ip 

Lh ^ 

llo 

0 

0 ) 


p+ 

H"! 


{h + -i) 


Ih J Ip 

1 -1 0 


( 68 ) 


[Lk]MJ] 

where is the parity function equal to one if n is even and zero if n is odd. We also use 
the following product of a nine-j and a three-j 


2 3p 

2 jh 

1 J' J 


'J 

0 0 0 


(_)Jp+h+2 

Vq [ip]MJ][J'] 

Xp + Xh + sJ + / jp jh J 


X 


V J + 


1 _1 0 

2 2 


(69) 


where Xp = + \)- After some work, we dually arrive at the matrix elements 

written in the form 


{p\\C?\\h) = ■^Aj{ph) drrR*{r)jj{qr)Rh{r) 


+ ^Bj{ph) drr‘^R*{r)jj{qr) 


dRhjr) 

dr 


^ ■? _ I fOO 

{p\\C?\\h) = — Y. Ajj,{ph)- / dr R*{r)jj,{qr)Rh{r) 

^ sA±i Q do ^ 

+ ^ Y drr^Rpir)jj,iqr)^^^^ 


(70) 


(71) 


where J' = J + s, s = ±1, and jji{z) is the derivative of the Bessel function. We have 
dedned the following coefdcients 


Aj{ph) 

— P + 

Bj{ph) 

= 

Ajj'iph) 

= P+ 

Bjj'iph) 

= P+ 

H"! 


[jyih{ih+i)bjj- Y 

s=±l rJ 


aj 






Mm] y 


jp jh J 


-d 0 


aj 


ro' 


hjJAi (xp+xk+< 

where the factors aj and bjj' are dedned as follows: 

2 y dp \ 

5jj' = \ \ ih jh 

1 J' J 


jp jh J 

1 _1 0 

r» r» W 


/ J' 1 

ih O ip 

1 -1 0 


(72) 

(73) 

(74) 

(75) 

(76) 

(77) 
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B PWIA 


We follow the approach of [27]. In PWIA the inclusive response functions can be written as 
an integral over the missing momentum p = — q, with p' = 


R 


PWIA 

K 




^ dpp f d(j)WKip\p)Mh{p) 

p' — q\ Jo 


ih q 

of the scalar momentum distribution for each occupied shell 

Mh{p) = -^^^^\Rh{p)\‘^ , 

where Rh{p) is the radial wave function in momentum space. The single¬ 
responses wk{p',p), for K = CC,CL, LL,T,T' are readily computed 
Eqs. (24,25), and axial-vector, Eqs. (28,29,30), current components 


(78) 


(79) 

nucleon exclusive 
using the vector. 


Wcc 


(80) 

V 

^cc 

= + 

(81) 

Wcc 

= Cow^ + C'^dl 

(82) 

WCL 

= + 

(83) 

WcL 

^ V 

= --Wcc 

(84) 

WcL 


(85) 

wll 

= wIl + w^l 

(86) 

wL 

= (^) wlc 

(87) 

wtp 

= 

(88) 

Wt 

= IV^ + 

(89) 


= 2^[\^+e,vi 

(90) 


= 2C(^ 

(91) 

Wt' 

= 

(92) 


Since these functions do not depend on the azimuthal angle of p, namely 0, the response 
functions in Eq. (78) are reduced to an integral over the missing momentum p, which is 
performed numerically. 


C The relativistic Fermi gas 

Here we summarize the expressions for (z/^,/i“) reactions in the REG model. We follow [37], 
where the expressions were written to leading order. Here we write the full results for the 
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non-Pauli blocked regime of interest in this work. The nuclear response functions are written 
as 

Rk = NAoUkIrfcW, K = CC, CL, LL, T, T', (93) 

where N is the neutron number, 

Ao = ^^3 . (94) 

rUNTIpK 

Here rjF = kp/mN, = \j^ + Vf ~ 1 - Eq. (93) fuFcii’) is the scaling function of the 
RFG 

fRFcW = ^(1 - ^^)^(1 - ^^) (95) 

and Ip is the scaling variable given in Eq. (55). Finally, we give the single-nucleon responses 
Uk- For K = CC we have 


where 


Ucc 


TtV 

^cc 


U, 


K 


V 

CC 
2 


r 


{2C^Y + 



(96) 

(97) 


A = -^^(1 - r ) 


kWI + ^ + y(l 


(98) 


and we have written the axial-vector response as the sum of conserved (c.) plus non conserved 
(n.c.) parts. 




- 

(99) 



- — . 

(100) 

For K = CL, LL we have 




UcL 

= U^,+ { 


(101) 

Ull 

= uL + ( 


(102) 

where the vector and conserved axial-vector parts are determined by current conservation as 


tC — 
'-'cl — 


(103) 


K). = 

-- (Go) 

K ^ /c. 

(104) 


- 

'-'ll — 

—U^ 

(105) 


LC. - 


(106) 
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while the n.c. parts are 




r 


/ 2 
A 



T 


/ 2 
A 


Finally the transverse responses are given by 


with 


[/^ = [/^ + U4 


= 

= 

Ut' = 


2T{2G\,f + + 

1 + r 

2{l +t)G\ + G\A 
2Ga{2G\j)Jt{1 + t) [1 + A] 


A = 


r 


1 + r 


2k 
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